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I It is proposed, that color-singlet gluon-clusters can be formed in hadrons as a con- 

t — 

I sequence of self-organized criticality (SOC) in systems of interacting soft gluons, and 

00 ' 

0^ I that the properties of such spatiotemporal complexities can be probed experimentally 

Q^l by examing inelastic diffractive scattering. Theoretical arguments and experimental 

p ,' 

Q ' evidences supporting the proposed picture are presented — together with the result of 

a systematic analysis of the existing data for inelastic diffractive scattering processes 
^ ' performed at different incident energies, and/or by using different beam-particles. It 



is shown that the size- and the lifetime-distributions of such gluon-clusters can be di- 
rectly extracted from the data, and the obtained results exhibit universal power-law 
behaviors — in accordance with the expected SOC-fingerprints. 
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1. Interacting soft gluons in the small-x^ region of DIS 



A number of striking phenomena have been observed in recent deep-inelastic electron-proton 
scattering (DIS) experiments in the small-x^ region. In particular it is seen, that the 
contribution of the gluons dominatesS, and that large-rapidity-gap (LRG) events existii. 
The latter shows that the virtual photons in such processes may encounter "colorless objects" 
originating from the proton. 

The existence of LRG events in these and otheri'i scattering processes have attracted 
much attention, and there has been much discussionHlSl on problems associated with the 
origin and/or the properties of such "colorless objects". Reactions in which "exchange" of 
such "colorless objects" dominate are known in the literatureii'i as "diffractive scattering 
processes". While the concepts and methods used by different authors in describing such 
processes are in general very much different from one another, all the authors (experimen- 
talists as well as theorists) seem to agree on the following (see also Refs. P-p!T[l): 
(a) Interacting soft gluons play a dominating role in understanding the phenomena in the 
small-xs region of DIS in general, and in describing the properties of LRG events in par- 
ticular, (b) Perturbative QCD should be, and can be, used to describe the LRG events 
associated with high transverse-momentum (p±) jets which have been observed at HERaI 
and at the Tevatroni. Such events are, however, rather rare. For the description of the bulk 
of LRG events, concepts and methods beyond the perturbative QCD, for example Pomeron 
ModelJl based on Regge Phenomenology, are needed. It has been suggested a long time ago 
(see the first two papers in Ref.i) that, in the QCD language, "Pomeron-exchange" can be 
interpreted as "exchange of two or more gluons" and that such results can be obtained by 
calculating the corresponding Feynman diagrams. It is generally felt that non-perturbative 
methods should be useful in understanding "the small-x phenomena", but the question, 
whether or how perturbative QCD plays a role in such non-perturbative approaches does 
not have an unique answer. 

In a recent LetteB, we proposed that the "colorte. objects" which play the dominating 
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role in LRG events are color-singlet gluon-clusters due to self-organized criticality, and that 
optical-geometrical concepts and methods are useful in examing the space-time properties 
of such objects. 

The proposed picture^ is based on the following observation: In a system of soft gluons 
whose interactions are not negligible, gluons can be emitted and/or absorbed at any time 
and everywhere in the system due to color-interactions between the members of the system 
as well as due to color-interactions of the members with gluons and/or quarks and antiquarks 
outside the system. In this connection it is important to keep in mind that gluons interact 
directly with gluons and that the number of gluons in a system is not a conserved quantity. 
Furthermore, since in systems of interacting soft-gluons the "running-coupling-constant" is 
in general greater than unity, non-perturbation methods are needed to describe the local 
interactions associated with such systems. That is, such sytems are in general extremly 
complicated, they are not only too complicated (at least for us) to take the details of local 
interactions into account (for example by describing the reaction mechanisms in terms of 
Feynman diagrams), but also too complicated to apply well-known concepts and methods 
in conventional Equilibrium Statistical Mechanics. In fact, the accumulated empirical facts 
about LRG events and the basic properties of gluons prescribed by the QCD are forcing us 
to accept the following picture for such systems: 

A system of interacting soft gluons can be, and should be considered as an open dynamical 
complex system with many degrees of freedom, which is in general far from equilibrium. 

In our search for an appropriate method to deal with such complex systems, we are led 
to the following questions: Do we see comparable complex systems in Nature? If yes, what 
are the characteristic features of such systems, and what can we learn by studying such 
systems? 
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2. Characteristic features of open dynamical complex systems 



Open, dynamical, complex systems which are in general far from equilibrium are not difficult 
to find in Nature — at least not in the macroscopic world! Such systems have been studied, 
and in particular the following have been observed by Bak, Tang and Wiesenfeld (BTW) 
some time ago0: This kind of complex systems may evolve to self-organized critical states 
which lead to fiuctuations extending over all length- and time-scales, and that such fiuc- 
tuations manifest themselves in form of spatial and temporal power-law scaling behaviors 
showing properties associated with fractal structure and fiicker noise respectively. 

To be more precise, and many other authorslll proposed, and demonstrated by 

numerical simulations, the following: Open dynamical complex systems of locally interacting 
objects which are in general far from equilibrium can evolve into self-organized structures of 
states which are barely stable. A local perturbation of a critical state may "propagate" , in 
the sense that it spreads to (some) nearest neighbors, and than to the next-nearest neighbors, 
and so on in a "domino effect" over all length scales, the size of such an "avalanche" can be 
as large as the entire system. Such a "domino effect" eventually terminates after a total time 
T, having reached a final amount of dissipative energy and having effected a total spatial 
extension S. The quantity S is called by BTW the "size", and the quantity T the "lifetime" 
of the avalanche — named by BTW a "cluster" (hereafter referred to as BTW-cluster or 
BTW-avalanche). As we shall see in more details later on, it is of considerable importance 
to note that a BTW-cluster cannot, and should not be identified with a cluster in the usual 
sense. It is an avalanche, not a static object with a fixed structure which remains unchanged 
until it decays after a time-interval (known as the lifetime in the usual sense). 

In fact, it has been shown0'Ei that the distribution (-Ds) of the "size" (which is a measure 
of the dissipative energy, S) and the distribution (-Dt) of the lifetime (T) of BTW-clusters 
in such open dynamical complex systems obey power-laws: 

DsiS) ~ (1) 
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DriT) ~ T-^ (2) 

where n and u are positive real constants. Such spatial and temporal power-law scaling 
behaviors can be, and have been, considered as the universal signals — the "fingerprints" — 
of the locally perturbed self-organized critical states in such systems. It is expectedliilii that 
the general concept of self-organized criticality (SOC), which is complementary to chaos, 
may be the underlying concept for temporal and spatial scaling in a wide class of open 
n on- equilibrium complex systems — although it is not yet known how the exponents of such 
power laws can be calculated analytically. 

SOC has been observed in a large number of open dynamical complex systems in non- 
equilibriumSilSii among which the following examples are of particular interest, because 
they illuminate several aspects of SOC which are relevant for the discussion in this paper. 

First, the well known Gutenberg- Richter law0'§ for earthquakes special case of 
Eq.(l): In this case, earthquakes are BTW-clusters due to SOC. Here, S stands for the 
released energy (the magnitude) of the observed earthquakes. Ds{S) is the number of 
earthquakes at which an energy S is released. Such a simple law is known to be valid for all 
earthquakes, large (up to 8 or 9 in Richter scale) or small! We note, the power-law behavior 
given by the Gutenberg-Richter law implies in particular the following. The question "How 
large is a typical earthquake?" does not make sense! 

Second, the sandpile experimentsllilii which show the simple regularities mentioned in 
Eqs.(l) and (2): In this example, we see how local perturbation can be caused by the 
addition of one grain of sand (note that we are dealing with an open system!). Here, we can 
also see how the propagation of perturbation in form of "domino effect" takes place, and 
develops into BTW-clusters/avalanches of all possible sizes and durations. The size- and 
duration-distributions are given by Eqs.(l) and (2) respectively. This example is indeed a 
very attractive one, not only because such experiments can be, and have been performed in 
laboratorieJii, but also because they can be readily simulated on a pdiilii. 

Furthermore, it has been pointed out, and demonstrated by simple mo that 
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the concept of SOC can also be applied to Biological Sciences. It is amazing to see how 
phenomena as complicated as Life and Evolution can be simulated by simple models such 
as the "Game of Life"0 and the "Evolution Model" ffl. 

Having seen that systems of interacting soft-gluons are open dynamical complex systems, 
and that a wide class of open systems with many degrees of freedom in the macroscopic world 
evolve to self-organized critical states which lead to fluctuations extending over all length- 
and time-scales, it seems natural to ask the following: Can such states and such fluctuations 
also exist in the microscopic world — on the level of quarks and gluons? In particular: Can 
SOC be the dynamical origin of color-singlet gluon-clusters which play the dominating role 
in inelastic diffractive scattering processes? 

3. SOC in inelastic diffractive scattering processes? 

Because of the special role played by "the colorless objects" in inelastic diffractive scattering, 
and the possible relations between such objects and color-singlet gluon-clusters which can 
be formed in systems of interacting soft gluons, it should be of considerable interest to study 
the questions mentioned at the end of the last section, as well as in the title of this section. A 
simple and effective way of answering them, is to check whether the characteristic properties 
of SOC, in particular the SOC- "fingerprints" mentioned in Eqs.(|l]) and (^) show up in the 
relevant experiments. In order to perform such a comparison, we need to extract the spatial 
and the temporal distributions of the gluon-clusters. 

What are such "colorless objects"? Is it possible that the colorless objects which are 
associated with the proton-target and which play the dominating role in inelastic diffractive 
scattering processes are BTW-clusters which exist due to SOC in systems of interacting soft 
gluons? Can we examine the properties of such colorless objects by studying the final states 
of the above-mentioned scattering processes? 

To answer these questions, it is useful to recall the following: As color- singlets, such 
colorless objects can exist inside and/or outside the proton, and the interactions between 



5 



such color-singlets as well as those between such objects and "the mother proton" should 
be of Van der Waals type. Hence it is expected that such a colorless object can be readily 
separated as an entire object from the mother proton in scattering processes in which the 
momentum-transfer is sufficient to overcome the binding energy due to the Van der Waals 
type of interactions. This means, in inelastic diffractive scattering the beam-particle (which 
is the virtual photon 7* in DIS) should have a chance to encounter on of the color-singlet 
gluon-clusters. For the reasons mentioned above, the struck colorless object can simply 
be "knocked out" and/or "carried away" by the beam-particle in such a collision event. 
Hence, it seems that the question whether "the colorless objects" are indeed BTW-clusters is 
something that can be answered experimentally. In this connection we recall that, according 
to the general theory of soci'i, the size of a BTW-cluster is characterized by its dissipative 
energy, and in case of systems of interacting soft gluons associated with the proton, the 
dissipative energy carried by the BTW-cluster should be proportional to the energy fraction 
(xp) carried by the colorless object. Hence, if the colorless object can indeed be considered 
as a BTW-cluster due to SOC, we should be able to obtain information about the size- 
distribution of such color-singlet gluon-clusters by examing the xp-distributions of LRG 
events in the small-x^ region of DIS. 

Having this in mind, we now take a closer look at the measured! "diffractive structure 
function" F^^^\(3,Q^;xp) = J dtF^^^\(3,Q^; xp,t). Here, we note that F^^^\f3,Q^; xp,t) 
IS relatedi^i to the differential cross-section for large-rapidity-gap events 

dpdQ^dxpdt = -W^'-y+Y^^' (^'^ 

in analogy to the relationship between the corresponding quantities [namely d^a /{dxp dQ^) 
and F2{xb, Q"^)] for normal deep-inelastic electron-proton scattering events 

''"'■i-y^'4)F.M''y (4) 



dxpdQ^ xbQ^ 2 
The kinematical variables, in particular (3, Q^, xp and xb (in both cases) are directly 
measurable quantities, the definitions of which are shown in Fig.l together with the corre- 
sponding diagrams of the scattering processes. We note that, although these variables are 
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Lorentz-invariants, it is sometimes convenient to interpret them in a "fast moving frame", 
for example the electron-proton center-of-mass frame where the proton's 3-momentum P is 
large (i.e. its magnitude \P\ and thus the energy = (|-Pp + M^)^/^ is much larger than 
the proton mass M). While characterizes the virtuality of the space-like photon 7*, xb 
can be interpreted, in such a "fast moving frame" (in the framework of the celebrated parton 
model), as the fraction of proton's energy P° (or longitudinal momentum \P\) carried by 
the struck charged constituent. 

We recall, in the framework of the parton model, F2{xb, Q'^)/xb for "normal events" can 
be interpreted as the sum of the probability densities for the above-mentioned 7* to interact 
with a charged constituent of the proton. In analogy to this, the quantity F-^^^\f3, Q'^; Xp)/(3 
for LRG events can be interpreted as the sum of the probability densities for 7* to interact 
with a charged constituent which carries a fraction (3 = xb/xp of the energy (or longitudinal 
momentum) of the colorless object, under the condition that the colorless object (which we 
associate with a system of interacting soft gluons) carries a fraction of proton's energy 
(or longitudinal momentum). We hereafter denote this charged-neutral and color- neutral 
gluon-system by Cq (in Regge pole mo delH this object is known as the "pomeron"). Hence, 
by comparing Eq. (3) with Eq. (4) and by comparing the two diagrams shown in Fig. 1(a) 
and Fig. 1(b), it is tempting to draw the following conclusions: 

The diffractive process is nothing else but a process in which the virtual photon 7* 
encounters a Cg, and P is nothing else but the Bj or ken- variable with respect to Cg (this is 
why it is called xbc in Ref . [p!0[| ) . This means, a diffractive e~p scattering event can be 
envisaged as an event in which the virtual photon 7* collides with "a Cg-target" instead of 
"the proton-target". Furthermore, since Cg is charge-neutral, and a photon can only directly 
interact with an object which has electric charges and/or magnetic moments, it is tempting 
to assign Cg an electro-magnetic structure function F2(/5, Q^), and study the interactions 
between the virtual photon and the quark(s) and antiquark(s) inside Cg. In such a picture 
(which should be formally the same as that of Regge pole mo delsi, if we would replace the 
Cg's by "pomerons") we are confronted with the following two questions: 
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First, is it possible and meaningful to discuss the xp-distributions of the Cq's without 
knowing the intrinsic properties, in particular the electromagnetic structures, of such ob- 
jects? 

Second, are gluon-clusters hadron-like, such that their electromagnetic structures can be 
studied in the same way as those for ordinary hadrons? 

Since we wish to begin the quantitative discussion with something familiar to most of the 
readers in this community, and we wish to differentiate between the conventional-approach 
and the SOC-approach, we would like to discuss the second question here, and leave the 
first question to the next section. We recall (see in particular the last two papers in Ref.^, 
in order to see whether the second question can be answered in the affirmative, we need to 
know whether F^^^\l3, Q^; xp) can be factorized in the form 



Here, fdxp) plays the role of a "kinematical factor" associated with the "target Cq", and 
Xp is the fraction of proton's energy (or longitudinal momentum) carried by Cq. [We could 
call fc{xp) "the Cg-flux" — in exactly the same manner as in Regge pole mo delJ, where it 
is called "the pomeron flux".] F2(/5, Q^) is "the electro-magnetic structure function of Cq" 
[the counterpart of F2{xb,Q'^) of the proton] which — in analogy to proton (or any other 
hadron) — can be expressed as 



where q'f(^f) stands for the probability density for 7* to interact with a quark (antiquark) 
of flavor i and electric charge which carries a fraction (3 of the energy (or longitudinal 
momentum) of Cq. It is clear that Eq.(6) should be valid for all xp- values in this kinematical 
region, that is, both the right- and the left-hand-side of Eq. (6) should be independent of the 
energy (momentum) carried by the "hadron" Cg. 

Hence, to find out experimentally whether the second question can be answered in the 
affirmative, we only need to check whether the data are in agreement with the assumption 



(5) 



P 



(6) 
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that F2(/3, Q^) prescribed by Eqs.(5) and (6) exists. For such a test, we take the existing 
data! and plot log[F2^*'^''(/9, Q^; xp)//5] against log/3 for different xp-values. We note, under 
the assumption that the factorization shown in Eq.(5) is valid, the /5-dependence for a 
given in such a plot should have exactly the same form as that in the corresponding 
log[F2(/3, Q^)//^] vs log/9 plot; and that the latter is the analog of log[F2(xB, vs 
logxp plot for normal events. In Fig.2 we show the result of such plots for three fixed Q^- 
values (3.5, 20 and 65 GeV^, as representatives of three different ranges in Q^). Our goal is 
to examine whether or how the /3-dependence of the function given in Eq.(6) changes with 
xp. In principle, if there were enough data points, we should, and we could, do such a plot 
for the data-sets associated with every xp-value. But, unfortunately there are not so much 
data at present. What we can do, however, is to consider the /^-distributions in different 
xp-bins, and to vary the bin-size of xp, so that we can explicitly see whether/how the shapes 
of the /^-distributions change. The results are shown in Fig.2. The /^-distribution in the first 
row, corresponds to the integrated value F^{/3, Q"^) shown in the literatureSi. Those in the 
second and in the third row are obtained by considering different bins and/or by varying 
the sizes of the bins. By joining the points associated with a given xp-interval in a plot 
for a given Q^, we obtain the /3-distribution for a Cg carrying approximately the amount 
of energy xpP^, encountered by a photon of virtuality Q^. Taken together with Eq.(6) we 
can then extract the distributions q^iP, Q^) and q^^P, Q^) for this Q^-value, provided that 
F2{(3,Q'^)/ (3 is independent of xp. But, as we can see in Fig.2, the existing datai'i show 
that the xp-dependence of this function is far from being negligible! Note in particular 
that according to Eq.(^), by choosing a suitable /p(xp) we can shift the curves for different 
Xp- values in the vertical direction (in this log-log plot); but we can never change the shapes 
of the j3- distributions which are different for different xp- values! 

In order to see, and to realize, the meaning of the xp-dependence of the distributions 
of the charged constituents of Cg expressed in terms of F2{P,Q'^)/ (3 in LRG events [see 
Eqs.(5) and (6)], let us, for a moment, consider normal deep-inelastic scattering events 
in the xp-region where quarks dominate (xp > 0.1, say). Here we can plot the data for 
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log[F2{x B , Q"^) / X b] as a function of logx^ obtained at different incident energies (P^'s) 
of the proton. Suppose we see, that at a given Q^, the data for x^-distributions taken 
at different values of are very much different. Would it still be possible to introduce 
F2{xb,Q'^) as "the electro-magnetic structure function" of the proton, from which we can 
extract the xs-distribution of the quarks qi{xBiQ'^) at a given Q^l The fact that it is not 
possible to assign an xp-independent structure function F2{P,Q^)/ P to Cq which stands for 
the "pomeron", and whose "flux" fc{xp) is expected to be independent of P and Q^, deserves 
to be taken seriously. It strongly suggest that the following picture cannot be true: "There 
exists a universal colorless object (call it pomeron or Cq or something else) the exchange 
of which describes diffractive scattering in general and DIS off proton in particular. This 
object is hadron-like in the sense that it has not only a typical size and a typical lifetime, 
but also a typical electromagnetic structure which can e.g. be measured and described by 
an "electromagnetic structure function". 

In summary of this section, we note that the empirical facts mentioned above show that 
no energy-independent electromagnetic strcture function can be assigned to the expected 
universal colorless object Cq. This piece of experimental fact is of considerable importance, 
because it is the first indication that, if there is a universal "colorless object", this object 
cannot be considered as an ordinary hadron. In other words, it has to be something else! 
In fact, as we shall see below, this property is closely related to the observation that such 
an object cannot have a typical size, or a typical lifetime. The final answer to the question 
mentioned in the title of this section will be presented in Section ^ 

4. Distributions of the gluon-clusters 

After having seen that the existing data does not allow us to assign an energy-independent 
electromagnetic structure function to "the colorless object" such that the universal colorless 
object (cq) can be treated as an ordinary hadron, let us now come back to the first question 
in Section |^, and try to find out whether it is never-the-less possible, and meaningful, to 
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talk about the xp-distribution of Cq. As we shall see in this section, the answer to this 
question is Yes! Furthermore, we shall also see, in order to answer this question in the 
affirmative, we do not need the factorization mentioned in Eq.(5), and we do not need to 
know whether the gluon-clusters are hadron-like. But, as we have already mentioned above, 
it is of considerable importance to discuss the second question so that we can understand 
the origin and the nature of the Cq's. 

In view of the fact that we do use the concept "distributions of gluons" in deep-inelastic 
lepton-hadron scattering, although the gluons do not directly interact with the virtual pho- 
tons, we shall try to introduce the notion "distribution of gluon-clusters" in a similar manner. 
In order to see what we should do for the introduction of such distributions, let us recall the 
following: 

For normal deep-inelastic e~p collision events, the structure function F2{xb,Q^) can be 
expressed in term of the distributions of partons, where the partons are not only quarks 
and antiquarks, but also gluons which can contribute to the structure function by quark- 
antiquark pair creation and annihilation. In fact, in order to satisfy energy-momentum- 
conservation (in the electron-proton system), the contribution of the gluons Xgg{xg, Q^) has 
to be taken into account in the energy-momentum sum rule for all measured Q^-values. Here, 
we denote by g{xg, Q^) the probability density for the virtual photon 7* (with virtuality Q^) 
to meet a gluon which carries the energy (momentum) fraction Xg of the proton, analogous 
to qi{xB,Q'^) [or qi{xB-,Q'^)] which stands for the probabihty density for this 7* to interact 
with a quark (or an antiquark) of flavor i and electric charge Cj which carries the energy 
(momentum) fraction xb oi the proton. We note, while both xb and Xg stand for energy 
(or longitudinal momentum) fractions carried by partons, the former can be, but the latter 
cannot be directly measured. 

Having these, in particular the energy-momentum sum rule in mind, we immediately see 
the following: In a given kinematical region in which the contributions of only one category 
of partons (for example quarks for Xb > 0.1 or gluons for Xb < 10^^) dominate, the structure 
function F2{xb, Q^) can approximately be related to the distributions of that particular kind 
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of partons in a very simply manner. In fact, the expressions below can be, and have been, 
interpreted as the probability-densities for the virtual photon 7* (with virtuality Q'^) to meet 
a quark or a gluon which carries the energy (momentum) fraction xb ot Xg respectively. 

F2{xb,Q'' 



Xb 



F2{xb,Q 



2^ 



Xg 



Y.^'i1ii^B,Q'^) or 



9{x,,Q'). (7) 



The relationship between qi{xB,Q^), g{xg,Q^) and F2{xb,Q^) as they stand in Eq.(^ are 
general and formal (this is the case especially for that between g and -F2) in the following 
sense: Both qi{xB, Q^) and g{xg, Q"^) contribute to the energy-momentum sum rule and both 
of them are in accordance with the assumption that partons of a given category (quarks or 
gluons) dominate a given kinematical region (here xb > 0.1 and xb < 10~^ respectively). 
But, neither the dynamics which leads to the observed Q^-dependence nor the relationship 
between Xg and xb are given. This means, without further theoretical inputs, the simple 
expression for g{xg,Q'^) as given by Eq.(7) is practically uselessl 

Having learned this, we now discuss what happens if we assume, in diffractive lepton- 
nucleon scattering, the colorless gluon-clusters (cg's) dominate the small-Xs region {xb < 
10"^, say). In this simple picture, we are assuming that the following is approximately true: 
The gluons in this region appear predominately in form of gluon clusters. The interaction 
between the struck Cg and the rest of the proton can be neglected during the 7-Cq collision 
such that we can apply impuls-approximation to the Cq's in this kinematical region. That 
is, here we can introduce — in the same manner as we do for other partons (see Eq.^), 
a probability density D5(xp|/5, Q^) for 7* in the diffractive scattering process to "meet" a 
Cq which carries the fraction xp of the proton's energy P*^ = (|-Pp -|- M'^Y^'^ ^ \P\ (where 
P is the momentum and M is the mass of the proton). In other words, in diffractive 
scattering events for processes in the kinematical region xp < 10^^, we should have, instead 
of g{xg,Q'^), the following: 

^f^M^ . DsiMO,Q-') . (8) 

Xp 
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Here, xpP^ is the energy carried by Cg, and (5 indicates the corresponding fraction carried by 
the struck charged constituent in Cq. In connection with the similarities and the differences 
between qi{xB, <5^), g{xB, Q"^) in and Ds{xp\l3, Q^) in (||), it is useful to note in particular 
the significant difference between Xg and xp, and thus that between the x^-distribution 
g{xg, Q"^) of the gluons and the xp-distribution Ds{xp\(3, Q"^) of the Cq's: Both Xg and xp are 
energy (longitudinal momentum) fractions of charge-neutral objects, with which 7* cannot 
directly interact. But, in contrast to Xg, xp can he directly measured in experiments, namely 
by making use of the kinematical relation 

and by measuring the quantities Q^, and W'^ in every collision event. Here, Q, and 
W stand respectively for the invariant momentum-transfer from the incident electron, the 
invariant-mass of the final hadronic state after the 7* — Cq collision, and the invariant mass 
of the entire hadronic system in the collision between 7* and the proton. Note that xp = 
(3xp, hence jS is also measurable. This means, in sharp contrast to g{xg,Q'^), experimental 
information on Ds{xp\l3,Q'^) in particular its xp-dependence can be obtained — without 
further theoretical inputs] 



5. The first SOC-fingerprint: Spatial scaling 

We mentioned at the beginning of Section ^, that in order to find out whether the concept 
of SOC indeed plays a role in diffractive DIS we need to check the fingerprints of SOC 
shown in Section |^, and that such tests can be made by examing the corresponding cluster- 
distributions obtained from experimental data. We are now ready to do this, because we have 
learned in Sections ^ and ^, that it is not only meaningful but also possible to extract x p- 
distributions from the measured diffractive structure functions, although the gluon-clusters 
cannot be treated as hadrons. In fact, as we can explicitly see in Eqs.(8) and (9), in order 
to extract the xp-dependence of the gluon-clusters from the data, detailed knowledge about 
the intrinsic structure of the clusters are not necessary. 
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Having these in mind, we now consider Ds as a function of xp for given values of (3 
and Q^, and plot F^^^\(3, Q^; xp)/xp against xp for different sets of (3 and Q^. The results 
of such log-log plots are shown in Fig. 3. As we can see, the datai suggest that the 
probability-density for the virtual photon 7* to meet a color-neutral and charged-neutral 
object Cq with energy (longitudinal momentum) fraction xp has a power-law behavior in 
Xp, and the exponent of this power-law depends very little on and p. This is to be 
compared with Ds{S) in Eq.( |l|), where S, the dissipative energy (the size of the BTW- 
cluster) corresponds to the energy of the system Cq. The latter is xpP^, where P° is the 
total energy of the proton. 

It means, the existing datai show that /^^(xpl/?, Q^) exhibits the same kind of power- 
law behavior as the size-distribution of BTW-clusters. This result is in accordance with 
the expectation that self-organized critical phenomena may exist in the colorless systems of 
interacting soft gluons in diffractive deep-inelastic electron-proton scattering processes. 

We note, up to now, we have only argued (in Section |1]) that such gluon-systems are 
open, dynamical, complex systems in which SOC may occur, and we have mentioned (in 
Section ^ the ubiquitousness of SOC in Nature. Having seen the experimental evidence 
that one of the "SOC-fingerprints" (which are necessary conditions for the existence of 
SOC) indeed exists, let us now take a second look at the colorless gluon-systems from a 
theoretical standpoint. Viewed from a "fast moving frame" which can for example be the 
electron-proton c.m.s. frame, such colorless systems of interacting soft gluons are part of 
the proton (although, as color- singlets, they can also be outside the confinement region). 
Soft gluons can be intermittently emitted or absorbed by gluons in such a system, as well 
as by gluons, quarks and antiquarks outside the system. The emission- and absorption- 
processes are due to local interactions prescribed by the well-known QCD-Lagrangian (here 
"the running coupling constants" are in general large, because the distances between the 
interacting colored objects cannot be considered as "short"; remember that the spatial 
dimension of a Cq can be much larger than that of a hadron!). In this connection, it is 
useful to keep the following in mind: Due to the complexity of the system, details about 
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the local interactions may be relatively unimportant, while general and/or global features 
— for example energy-flow between different parts (neighbors and neighbor's neighbors . . .) 
of the system — may play an important role. 

How far can one go in neglecting dynamical details when one deals with such open com- 
plex systems? In order to see this, let us recall how Bak and Sneppen0 succeeded in modeling 
some of the essential aspects of The Evolution in Nature. They consider the "fitness" of 
different "species", related to one another through a "food chain", and assumed that the 
species with the lowest fitness is most likely to disappear or mutate at the next time-step 
in their computer simulations. The crucial step in their simulations that drives evolution 
is the adaption of the individual species to its present environment (neighborhood) through 
mutation and selection of a fitter variant. Other interacting species form part of the envi- 
ronment. This means, the neighbors will be influenced by every time-step. The result these 
authors obtained strongly suggests that the process of evolution is a self-organized critical 
phenomenon. One of the essential simplifications they made in their evolution mo IS 
the following: Instead of the explicit connection between the fitness and the configuration of 
the genetic codes, they use random numbers for the fitness of the species. Furthermore, as 
they have pointed out in their papers, they could in principle have chosen to model evolution 
on a less coarse-grained scale by considering mutations at the individual level rather than 
on the level of species, but that would make the computation prohibitively difficult. 

Having these in mind, we are naturally led to the questions: Can we consider the creation 
and annihilation processes of colorless systems of interacting soft gluons associated with a 
proton as "evolution" in a microscopic world? Before we try to build models for a quan- 
titative description of the data, can we simply apply the existing evolution mo del^llli to 
such open, dynamical, complex systems of interacting soft-gluons, and check whether some 
of the essential features of such systems can be reproduced? 

To answer these questions, we now report on the result of our first trial in this direction: 
Based on the fact that we know very little about the detailed reaction mechanisms in such 
gluon-systems and practically nothing about their structures, we simply ignore them, and 
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assume that they are self-similar in space (this means, colorless gluon-clusters can be consid- 
ered as clusters of colorless gluon-clusters and so on). Next, we divide them in an arbitrary 
given number of subsystems Si (which may or may not have the same size) . Such a system is 
open, in the sense that neither its energy Si, nor its gluon-number has a fixed value. Since 
we do not know, in particular, how large the e^'s are, we use random numbers. As far the nj's 
are concerned, since we do not know how these numbers are associated with the energies in 
the subsystems Sj, except that they are not conserved quantities, we just ignore them, and 
consider only the e^'s. As in Ref. |]17| or in Ref.[jl8|, the random number of this subsystem 



as well as those of the fixedlll or random (see the first paper of Ref.[|l^]) neighbors will be 
changed at every time-step. Note, this is how we simulate the processes of energy flow due to 
exchange of gluons between the subsystems, as well as those with gluons/quarks/antiquarks 
outside the system. In other words, in the spirit of Bak and Sneppen0 we neglecting the 
dynamical details totally. Having in mind that, in such systems, the gluons as well as the 
subsystems (sj's) of gluons are virtual (space-like), we can ask: "How long can such a color- 
less subsystem Sj of interacting soft gluons exist, which carries energy ^j?" According to the 
uncertainty principle, the answer should be: "The time interval in which the subsystem Sj 
can exist is proportional to l/Si, and this quantity can be considered as the lifetime of Sj." 
In this sense, the subsystems of colorless gluons are expected to have larger probabilities 
to mutate because they are associated with higher energies, and thus shorter "lifetimes". 
Note that the basic local interaction in this self-organized evolution process is the emis- 
sion (or absorption) of gluons by gluons prescribed by the QCD-Lagrangian — although 
the detailed mechanisms (which can in principle be explicitly written down by using the 
QCD-Lagrangian) do not play a significant role. 

In terms of the evolution model0lli we now call Si the "species" and identify the corre- 
sponding lifetime Tj as the "fitness of Sj" . Because of the one-to-one correspondence between 
Ti and Ei, where the latter is a random number, we can also directly assign random numbers 
to the Tj's instead. From now we can adopt the evolution mo del0@ and note that, at the 
start of such a process (a simulation), the fitness on average grow, because the least fit are 
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always eliminated. Eventually the fitness do not grow any further on average. All gluons 
have a fitness above some threshold. At the next step, the least fit species (i.e. the most 
energetic subsystem Si of interacting soft gluons), which would be right at the threshold, 
will be "replaced" and starts an avalanche (or punctuation of mutation events), which is 
causally connected with this triggering "replacement". After a while, the avalanche will 
stop, when all the fitnesses again will be over that threshold. In this sense, the evolution 



goes on, and on, and on. As in Refs.[|l3l and [|1^], we can monitor the duration of every 
avalanche, that is the total number of mutation events in everyone of them, and count how 
many avalanches of each size are observed. The avalanches mentioned here are special cases 
of those discussed in Section |[ Their size- and lifetime-distributions are given by Eq.(l) 
and Eq.(2) respectively. Note in particular that the avalanches in the Bak-Sneppen model 
correspond to sets of subsystems Sj, the energies (cj) of which are too high "to be fit for the 
colorless systems of low-energy gluons". It means, in the proposed picture, what the virtual 
photon in deep-inelastic electron-proton scattering "meet" are those "less fit" one — those 
who carry "too much" energy. In a geometrical picture this means, it is more probable for 
such "relatively energetic" colorless gluons-clusters to be spatially further away from the 
(confinement region of) the proton. 

There exists, in the mean time, already several versions of evolution mo delsEi@ based on 
the original idea of Bak and Sneppen0 Although SOC phenomena have been observed in all 
these cases the slopes of the power-law distributions for the avalanches are different 

in different models — depending on the rules applied to the mutations. The values range 
from approximately —1 to approximately —2. Furthermore, these modelsElSlli seem to 
show that neither the size nor the dimension of the system used for the computer simulation 
plays a significant role. 

Hence, if we identify the colorless charge-neutral object Cq encountered by the virtual 
photon 7* with such an avalanche, we are identifying the lifetime of Cg with T, and the "size" 
(that is the total amount of dissipative energy in this "avalanche") with the total amount 
of energy of Cg. Note that the latter is nothing else but xpP^, where P° is the total energy 
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of the proton. This is how and why the S'-distribution in Eq. ([|) and the xp-distribution of 
Ds{xp\P,Q^) in Eq.(||) are related to each other. 

6. The second fingerprint: Temporal scaling 

In this section we discuss in more detail the effects associated with the time-degree-of- 
freedom. In this connection, some of the concepts and methods related to the two questions 
raised in Section ^ are of great interest. In particular, one may wish to know why the 
parton-picture does not work equally well for hadrons and for gluon-clusters. The answer 
is very simple: The time-degree of freedom cannot be ignored when we apply impulse- 
approximation, and the applicability of the latter is the basis of the parton-model. We 
recall that, when we apply the parton model to stable hadrons, the quarks, antiquarks and 
gluons are considered as free and stable objects, while the virtual photon 7* is associated 
with a given interaction-time Tint(<5^, a^s) characterized by the values and xb of such 
scattering processes. We note however that, this is possible only when the interaction- 
time Tint is much shorter than the corresponding time-scales related to hadron-structure (in 
particular the average propagation-time of color-interactions in hadron). Having these in 
mind, we see that, we are confronted with the following questions when we deal with gluon- 
clusters which have finite lifetimes: Can we consider the Cq's as "/ree" and "sto6/e" particles 
when their lifetimes are shorter than the interaction-time Tint(<5^, x^)? Can we say that a 
7*-Co collision process takes place, in which the incident 7* is absorbed by one a or a system 
of the charged constituents of Cg, when the lifetime T of Cq is shorter than Tint(<5^, a;^)? 

Since the notion "stable objects" or "unstable objects" depends on the scale which is 
used in the measurement, the question whether a Cq can be considered as a parton (in the 
sense that it can be considered as a free "stable object" during the 7*-Cq interaction) depends 
very much on on the interaction-time TintiQ"^ , x b) ■ Here, for given values of Q^, x^, and thus 
TintiQ"^ , X b) , only those Cg's whose lifetime (T's) are greater than TintiQ"^ ,xb) can absorb 
the corresponding 7*. That is to say, when we consider diffractive electron-proton scattering 
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in kinematical regions in which Cq's dominate, we must keep in mind that the measured 
cross-sections (and thus the diffractive structure function F^^^^) only include contributions 
from collision-events in which the condition T > Tmt{Q^,XB) is satisfied! 

We note that Tint can be estimated by making use of the uncertainty principle. In fact, 
by calculating in the above-mentioned reference frame, we obtain 
4|P| xb 



which implies that, for given \P\ and values. 



(10) 



Tint oc Xb, for Xb < 1- (11) 

This means, for diffractive e~p scattering events in the small-x^ region at given \P\ and 
values, xb is directly proportional to the interaction time Tint- Taken together with the 
relationship between Tint and the minimum lifetime T(mim) of the Cg's mentioned above, 
we reach the following conclusion: The distribution of this minimum value, T(min) of the 
Cq's which dominate the small-x^ {xb < 10^^, say) region can be obtained by examining 
the x^-dependence of F<^^^\f3,Q'^;xp)/(3 discussed in Eqs. (5), (6) and in Fig. 2. This is 
because, due to the fact that this function is proportional to the quark (antiquark) distribu- 
tions Qiiqf) which can be directly probed by the incident virtual photon 7*, by measuring 
F^^^\l3, Q^, xp)//3 as a function of xb = /3xp, we are in fact asking the following questions: 
Do the distributions of the charged constituents of Cq depend on the interaction time Tint, 
and thus on the minimum lifetime T(min) of the to be detected gluon-clusters ? We use 
the identity xb = I3xp and plot the quantity F.^^^\(3,Q'^; xp)/ (3 against the variable xp for 
fixed values of j3 and Q^. The result of such a log- log plot is given in Fig.4. It shows not 
only how the dependence on the time-degree-of-freedom can be extracted from the existing 
datai, but also that, for all the measured values of f3 and Q^, the quantity 

is approximately independent of /3, and independent of Q^. This strongly suggests that the 
quantity given in Eq.(0) is associated with some g'/o6a/ features of Cq — consistent with the 
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observation made in Section ^ which shows that it cannot be used to describe the structure 
of Cq. This piece of empirical fact can be expressed by setting ~ p{xb)- By 

taking a closer look at this log-log plot, as well as the corresponding plots for different sets 
of fixed (3- and Q^- values (such plots are not shown here, they are similar to those in Fig. 3), 
we see that they are straight lines indicating that p{xb) obeys a power-law. What does this 
piece of experimental fact tell us? What can we learn from the distribution of the lower 
limit of the lifetimes (of the gluon-systems Cq's)? 

In order to answer these questions, let us, for a moment, assume that we know the 
lifetime-distribution Dt{T) of the Cq's. In such a case, we can readily evaluate the integral 



and thus obtain the number density of all those clusters which live longer than the interaction 
time Ti^ti^B)- Hence, under the statistical assumption that the chance for a 7* to be absorbed 
by one of those Cg's of lifetime T is proportional to Dt(T) (provided that rint(Q^,a;B) < T, 
otherwise this chance is zero), we can then interpret the integral in Eq.(13) as follows: 
-^[Tint(Q^) xb)] oc p{xb) is the probability density for 7* [associated with the interaction-time 
TintixB)] to be absorbed by Cg's. Hence, 

Dt{xb) oc -^p{xb). (14) 



This means in particular, the fact that p{xb) obeys a power-law in implies that Dt{T) 
obeys a power-law in T. Such a behavior is similar to that shown in Eq.( H). In order to see 
the quality of this power- law behavior of Dt, and the quality of its independence of and 
f3, we compare the above-mentioned behavior with the existing datai. In Fig. 5, we show 
the log-log plots of d/dxB[p{xB)] against xb- We note that d/dxB[p{xB)] is approximately 
F^^^\(3, Q^; Xb/ I3)/{I3xb)- The quality of the power-law behavior of Dj- is explicitly shown 
in Fig. 5. 




(13) 



dXB 
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7. Q^-dependent exponents in the power-laws? 

We have seen, in Sections ^ and ^ that in diffractive deep-inelastic electron-proton scatter- 
ing, the size- and the lifetime-distributions of the gluon-clusters obey power-laws, and that 
the exponents depend very little on the variables (3 and Q^. We interpreted the power-law 
behaviors as the fingerprints of SOC in the formation processes of such clusters in form 
of BTW-avalanches. Can such approximately independence (or weak dependence) of the 
exponents on and (3 be understood in a physical picture based on SOC? In particular, 
what do we expect to see in photoproduction processes where the associated value for is 
zero? 

In order to answer these questions, let us recall the space-time aspects of the collision 
processes which are closely related to the above-mentioned power-law behaviors. Viewed 
in a fast moving frame (e.g. the c.m.s. of the colliding electron and proton), the states 
of the interacting soft gluons originating from the proton are self-organized. The colorless 
gluon-clusters caused by local perturbations and developed through "domino effects" are 
BTW-avalanches (see Sections |l| and ^, the size-distribution of which [see Eqs.(8) and (1)] 
are given by Fig.3. This explicitly shows that there are gluon-clusters of all sizes, because 
a power-law size-distribution implies that there is no scale in size. Recall that, since such 
clusters are color-singlets, their spatial extensions can be much larger than that of the proton, 
and thus they can be "seen" also outside the proton by a virtual photon originating from 
the electron. In other words, what the virtual photon encounters is a cloud of colorless 
gluon-clusters everyone of which is in general partly inside and partly outside the proton. 

The virtual photon, when it encounters a colorless gluon-cluster, will be absorbed by 
the charged constituents (quarks and antiquarks due to fluctuation of the gluons) of the 
gluon-system. Here it is useful to recall that in such a space-time picture, is inversely 
proportional to the transverse size, and is a measure of the interaction time [See Eqs. 
(10) and (11) in Section ^ of the virtual photon. It is conceivable, that the values for the 
cross-sections for virtual photons (associated with a given Q"^ and a given xb) to collide with 
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gluon-clusters (of a given size and a given lifetime) may depend on these variables. But, since 
the processes of self-organization (which produce such gluon-clusters) take place independent 
of the virtual photon (which originates from the incident electron and enters "the cloud" to 
look for suitable partners), the power-law behaviors of the size- and lifetime-distributions 
of the gluon-clusters are expected to be independent of the properties associated with the 
virtual photon. This means, by using 7*'s associated with different values of to detect 
clusters of various sizes, we are moving up or down on the straight lines in the log-log plots 
for the size- and lifetime distributions, the slopes of which do not change. In other words, the 
approximative Q^-independence of the slope is a natural consequence of the SOC picture. 

As far as the /3-dependence is concerned, we recall the results obtained in Sections |^ 
and ^, which explicitly show the following: The gluon-clusters (cq's) cannot be considered 
as hadrons. In particular, it is neither possible nor meaningful to talk about "the electro- 
magnetic structure of the gluon-cluster" . This suggests, by studying the /^-dependence of 
the "diffractive structure functions" we cannot expect to gain further information about the 
structure of the gluon-clusters or further insight about the reaction mechanisms. 

Having seen these, we try to look for measurable quantities in which the integrations over 
/? have already been carried out. A suitable candidate for this purpose is the differential 
cross-section 



Together with Eqs.(3) and (8), we see that this cross-section is nothing else but the effective 
/3-weighted xp-distribution Dsi^xplQ"^ , [3) of the gluon-clusters. Note that the weighting 
factors shown on the right- hand- side of Eq.(15) are simply results of QED! Next, we use the 
data! for ^2^*^^'' which are available at present, to do a log-log plot for the integrand of the 
expression in Eq.(15) as a function of xp for different values of (3 and Q^. This is shown 




(15) 
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in in Fig. 6a. Since the absolute values of this quantity depend very much, but the slope of 

the curves very little on j3, we carry out the integration as follows: We first fit every set 

of the data separately. Having obtained the slopes and the intersection points, we use the 

obtained fits to perform the integration over /3. The results are shown in the 

/ 1 rfV^ \ 

— -77^ — versus log (xp) 
\xp dQ'^ axp J 

plots of Fig. 6b. These results show the Q^-dependence of the slopes is practically negligible, 
and that the slope is approximately —1.95 for all values of Q^. 

Furthermore, in order to see whether the quantity introduced in Eq.(15) is indeed useful, 
and in order to perform a decisive test of the Q^-independence of the slope in the power-law 
behavior of the above-mentioned size-distributions, we now compare the results in deep- 
inelastic scattering! with those obtained in photoproductionEl, where LRG events have also 
be observed. This means, as in diffractive deep-inelastic scattering, we again associate the 
observed effects with colorless objects which are interpreted as system of interacting soft 
gluons originating from the proton. In order to find out whether it is the same kind of 
gluon-clusters as in deep-inelastic scattering, and whether they "look" very much different 
when we probe them with real (Q^ = 0) photons, we replot the existing da / dM^ dataS for 
photoproduction experiments performed at different total energies, and note the kinematical 
relationship between M^, W"^ and xp (for Q"^ < and |t| < Mj): 

- ^ (16) 
The result of the corresponding 

log ( ; — I versus \og(xp) 

\xp dxp J 

plot is shown in Fig. 7. The slope obtained from a least-square fit to the existing data0 is 
-1.98 ±0.07. 

The results obtained in diffractive deep-inelastic electron-proton scattering and that for 
diffractive photoproduction strongly suggest the following: The formation processes of gluon- 
clusters in the proton is due to self-organized criticality, and thus the spatial distributions 
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of such clusters — represented by the xp-distribution — obey power-laws. The exponents 
of such power-laws are independent of . Since can be interpreted as a measure 

for the transverse size of the incident virtual photon, the observed Q^-independence of the 
exponents can be considered as further evidence for SOC — in the sense that the self- 
organized gluon-cluster formation processes take place independent of the photon which is 
"sent in" to detect the clusters. 

Having these results, and the close relationship between real photon and hadron in mind, 
we are immediately led to the following questions: What shall we see, when we replace the 
(virtual or real) photon by a hadron — a proton or an antiproton? (See in this connection 
Fig. 8, for the notations and the kinematical relations for the description of such scattering 
processes.) Should we not see similar behaviors, if SOC in gluon-systems is indeed the 
reason for the occurrence of colorless gluon-clusters which can be probed experimentally in 
inelastic diffractive scattering processes? To answer these questions, we took a closer look at 
the available single diffractive proton-proton and proton-antiproton scattering data@@; and 
in order to make quantitative comparisons, we plot the quantities which correspond to those 
shown in Fig.^ and Fig.[7[ These plots are shown in Fig.|^ and Fig.p3. In Fig.|^a, we see 
the double differential cross-section {1/ xp)d?a / {dtdxp) at four different t-values. In Fig.^, 
we see the integrated differential cross-section {l/xp)da/dxp. Note that, here 

xp ^ Ml/s, (17) 

where ^/s is the total c.m.s. energy of the colliding proton-proton or antiproton-proton 
system. Here, the integrations of the double differential cross-section over t are in the 
ranges in which the corresponding experiments have been performed. (The extremely weak 
energy-dependence has been ignored in the integration.) The dashed lines in all the plots 
in Figs.|^ and ^3 stand for the slope —1.97 which is the average of the slope obtained from 
the plots shown in Figs.^ and |^. This means, the result shows exactly what we expect 
to see: The fingerprints of SOC can be clearly seen also in proton- and antiproton-induced 
inelastic diffractive scattering processes, showing that such characteristic features are indeed 
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universal and robust! 

We are thus led to the following conclusions. Color-singlet gluon-clusters can be formed 
in hadrons as a consequence of self-organized criticality (SOC) in systems of interacting 
soft gluons. In other words, "the colorless objects" which dominate the inelastic diffr active 
scattering processes are BTW-avalanches (BTW-clusters). Such color-singlet gluon-clusters 
are in general distributed partly inside and partly outside the confinement region of the 
"mother-hadron" . Since the interactions between the color-singlet gluon-clusters and other 
color singlet objects (including the target proton) should be of Van der Waals type, it is 
expected that such an object can be readily driven out of the above-mentioned confinement 
region by the beam-particle in geometrically more peripheral collisions. This is why we 
examined inelastic single-diffractive scattering processes at high energies in which virtual 
photon, real photon, proton, and antiproton are used as beam particles. This is also why we 
can extract the universal distributions of such color-singlet gluon-clusters directly from the 
data. In particular, the fact that is the energy fraction carried by the struck colorless 
gluon- cluster, and the fact that the xp-distributions are universal, it is tempting to regard 
such xp-distributions as "the parton-distributions" for diffractive scattering processes. Can 
we make use of such "parton-distributions" to describe and/or to predict the measurable 
cross-sections in inelastic diffractive scattering processes? This and other related questions 
will be discussed in Part II of the present paper. 
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FIGURES 

Fig. 1. The well-known Feynman diagrams (a) for diffractive and (b) for normal deep-inelastic 
electron-proton scattering are shown together with the relevant kinematical variables which describe 
such processes. 

Fig. 2. F2^^\p,Q^;xp)/l3 is plotted as a function of /? for given xp-intervals and for fixed 
Q^-values. The data are taken from Ref.[^. The lines are only to guide the eye. 

Fig. 3. F^^^\p,Q^;xp)/xp is plotted as a function of for different values of /? and Q^. The 
data are taken from Ref.|p. 

Fig. 4. F^^^\p,Q'^;xp)/P is plotted as a function oi xp in the indicated f3- and Q^-ranges. 
The data are taken from Ref.Q. 

Fig. 5. F^^^\p, Q^; xp / (5) / {I5xp) is plotted as a function of xp for fixed [3- and Q^-values. The 
data are taken from Ref.|^. 

Fig. 6. (a) {l/xp)(fa^/d(3dQ'^dxp in units of GeV"^ is plotted function of xp in different 
bins of (3 and . The data are taken from Ref.@. (b) {l/xp)d'^a^ /dQ^dxp in units of GeV~^ is 
plotted as a function of xp in different bins of Q^. The data are taken from Ref.[p. 

Fig. 7. (l/xp)d(T/dxp for photoproduction ^ + p ^ X + p is plotted as a function of xp. The 
data are taken from Ref.[^. Note that the data in the second paper are given in terms of relative 
cross sections. Note also that the slopes of the straight-lines are the same. The two dashed lines 
indicate the lower and the upper limits of the results obtained by multiplying the lower solid line 
by cr^Qi = 154 it 16(stat.) it 32(syst.)^6. This value is taken from the third paper in Ref.|19]. 



Fig. 8. Diagrams for different single diffractive reactions, together with the definitions of the 
relevant kinematic variables. 
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Fig. 9. a) {l/xp)(Pa/dxpdt for single diffractive p + p ^ p + X and p + p ^ p + X reactions 
is plotted as a function of at different values of t and ^/s. The data are taken from Ref.m|. 
b) The integrated (with respect to two different |t|-ranges) differential cross section [1/ xp)da / dxp 
for single diffractive p + p ^ p + X and p + p ^ p + X reactions is plotted as a function of xp . 
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